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We introduce and develop an approach to realizing a topological phase transition and non-Abelian 
statistics with dynamically induced Floquet Majorana fermions (FMFs). When the periodic driv- 
ing potential does not break fermion parity conservation, FMFs can encode quantum information. 
Quasi-energy analysis shows that a stable FMF zero mode and two other satellite modes exist in 
a wide parameter space with large quasi-energy gaps, which prevents transitions to other Floquet 
states under adiabatic driving. We also show that in the asymptotic limit FMFs preserve non- 
Abelian statistics and, thus, behave like their equilibrium counterparts. 
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Introduction — Proposals of solid state [H-Q and cold 
atomic [§-flo| systems hosting Majorana fermions (MFs) 
have been a recent focus of attention. These systems 
present novel prospects for quantum computation since 
a widely separated pair of MF bound states, that for- 
mally correspond to zero-energy states of an effective 
Bogoliubov-de Genncs (BdG) Hamiltonian, forms a non- 
local fermionic state that is immune to local sources of 
decoherence. Moreover, MFs obey non-Abelian statistics 
and thus have potential for topological quantum informa- 
tion processing. Among the key si gna tures of MFs are 
a zero- bias resonance in tunneling jlll [l2j , half- integer 
conductance quantization [l3|, EH , and a 47r Josephson 
effect [THJl . Some of these predictions have already re- 
ceived possible experimental support (l6l - fl8j . 

Topological states of matter can be induced dynami- 
cally by time-periodic driving, the so-called Floquet ap- 
proach num. This brought to the agenda the new con- 
cept of Floquet Majorana fermion (FMFs) [Io[ • It turns 
out that even if the system is initially in the topologically 
trivial state, its Floquet version may exhibit topological 
properties. A realization of such states where they can 
be readily manipulated and precisely tuned in a wide pa- 
rameter space is therefore highly desirable. The natural 
questions for FMF systems are: whether they are robust 
and tunable, whether they can encode quantum informa- 
tion, and whether they follow non-Abelian statistics as 
for their equilibrium counterparts. Our study aims to 
answer these questions. 

We consider a generic platform to investigate non- 
Abelian statistics and potentially to realize topolog- 
ical quantum computation based on FMFs. The 
model is broadly applicable to both semiconductor- 
superconductors heterostructures with strongspin-orbit 
interaction and in-plane magnetic field 0, 13] , and to 
cold atomic systems where superconducting order is con- 
trolled by Feshbach resonances while spin-orbit coupling 
and Zeeman field effects are induced by an optical Ra- 
man transition [Ioj |. The latter realization is practically 
more promising since it allows a greater degree of con- 



trol. Furthermore, atomic condensates can be isolated 
thus suppressing dissipation on long time scales. 

We show, first, that if FMFs exist, they will exist at 
any instantaneous time. Therefore, FMFs can encode 
quantum information if the driving potential does not 
break fermion parity conservation. We study the quasi- 
energy spectrum of the problem analytically by using a 
rotating frame analysis in the limit that the frequency is 
large compared to the band width. We also perform ex- 
act numerical calculations which capture certain features 
of the spectrum beyond the rotating wave approximation 
(RWA). A broad range of parameters supporting FMFs is 
identified as a function of driving frequency uj and ampli- 
tude K for three specific driving scenarios: periodic mod- 
ulation of the chemical potential, or the Zeeman field. 
Finally, by using a two-time formalism 0, Hj], we show 
that FMFs follow the same non-Abelian statistics as their 
stationary counterparts. This conclusion stems from the 
observation that a generalized Floquet Berry matrix does 
not affect the non-Abelian statistics of FMFs since large 
quasi-energy gap ensures no transitions to other Floquet 
quasi-energy states in the adiabatic movement. 

Floquet Theorem for Majorana Fermion — Let us con- 
sider Floquet theory {24|. Suppose that the Hamilto- 
nian has an explicit time dependence H (t) = H(t + T) 
with period T = 2tt/uj, where w is the driving fre- 
quency. The solution of the Schrodinger equation can 
be described by a complete set of time-dependent state 



l*a(<)> 



f '\(f> a (t)}, where quasi-energies e a sat- 



isfy the equation [H(t) — idt]\4> a (t)) = e a \4> a (t)) and 
I = \<t>a(t+T)) are Floquet states (hereafter h = 1). 
The evolution operator U(t) = Tcxp(— i H(t)dt) has 
the following property 



U(t + T,t)\</> a (t)) =e 



(1) 



One can define an effective Hamiltonian H e s(t) through 
the relation [H[21| 



U(t + T,t) 



= e -iHeff «T 



(2) 



2 



with H e ff(t)\cf) a (t)) = e a \<p a (t)). We treat t as just a pa- 
rameter. The effective Floquet Hamiltonian is defined at 
each instantaneous time, and the topological properties 



21J. 



of each of these Hamiltonians is the same [2C 

If the system is described by a BdG Hamiltonian, the 
quasi-particle excitation spectrum will possess a particle- 
hole symmetry even if the time-dependent potential is 
added For any quasi-energy state \<f> e (t)) = 7 e (i)|0), 
the relation %(t) = 7_ e (i) is guaranteed. So, the zero 
quasi-energy state reveals the existence of a Floquet MF 
jlOI | . The full wavefunction for eo = can be written as 
|*o(*)> = e-^\Mt)} = \M*)) = 7o(*)|0>, with 70 (t) = 
7q (t) . Since quasi-energy is only defined in a interval of lo 
(e.g. from — ui/2 to w/2), another type of Floquet MF ex- 



ists at e = ±lu/2 with e 



-iwt/2 



lu/2 



0- 



From the argument above, we can show that if the zero 
quasi-energy state exists, a zero energy Floquet MF mode 
"f(t) exists at any instantaneous time t. The MF operator 
evolves in time periodically 7(t) = y(t + T); in general, it 
is different at different instantaneous times, j(t) ^ jit'). 

Quasi-Energy Spectrum and Floquet Majorana 
Fermion — To demonstrate the existence of FMFs 
consider a one dimensional wire with Rashba spin-orbit 
interaction Ago, Zceman splitting V z , and proximity- 
induced superconducting term A. The system can be 
described by a tight-binding Hamiltonian [fl @, EH : 



H a = 



i i 

+ A so ^ ( 6 l+it% _ 4+i| Cit + h - c ) ' ( 3 ) 



Here, i and a =J \^ denote fermion site and spin indices 
while Cjo-(ct-) are corresponding operators, r\ is the hop- 
ping term along the chain which yields a band width 
D = 4r), and [i is the chemical potential of the lat- 
tice model which is set to the particlc-hole symmetric 



point |25(. Note that Hamiltonian Eq. ^ is equally 



generic for a system of cold atoms ■ 

To add time dependence, it is natural to consider mod- 
ulating one of the parameters in Hq: the chemical poten- 
tial, the spin-orbit coupling [28|, or the Zeeman field. We 
first consider periodic modulation of the chemical poten- 
tial; the Hamiltonian is H(t) = Hq + H^(t) with 



Hfi(t) = Kcos(ujt) y^jh^ + n iX ), 



(4) 



where hi a — c ia Ci a . To calculate the quasi-energy, one 
can choose a basis in the rotating frame [l[ 

\{n ia }; m) = e - liL ^^^+n^^t > (5) 

where is the basis of the unperturbed system, 

and m labels the photon sector of the Floquet basis. The 



quasi-energy can be obtained by diagonalizing the Flo- 
quet operator H(t) — idt in this basis. The orthonor- 
mality condition of the Floquet states is only defined 
in an extended Hilbert space [H, so the inner product 
must include an extra time integral over a full period: 
((•|-)) = (l/T) j Q T dt(-\-). The matrix elements read 

(({n w };m\H(t) - id t \{n' ta };m')) 



T 



xe 



■ Ei("iT+™u)|{ n ' y\ e -i(m-m')ut 



(6) 



Since different photon sectors are separated by an energy 
gap of order u>, in the limit w>D, the admixture of pho- 
ton sectors can be neglected; this is in essence the rotat- 
ing wave approximation. Then, we can consider only the 
zero photon sector and obtain an effective Floquet Hamil- 
tonian by computing the m = m! = matrix element. 
The key point to notice is that only the superconducting 
term in ([3]) fails to commute with the chemical potential 
operator X^("-s ; t + ^ti)- Evaluation of Eq. within the 
RWA yields an effective Floquet Hamiltonian with ex- 
actly the same form as Hq with the pairing coupling A 
effectively renormalized to 



A eff = AJ (2K/cj). 



(7) 



(Jo(x) is the zero order Bessel function of the first kind.) 

We conclude from Eq. ([7]) that in Floquet systems one 
may induce a topological phase transition dynamically. 
Indeed, recall that the regime for a topological supercon- 
ducting phase of Hq, which supports MFs, requires the 
condition V, 2 > A 2 + (fx + 2rf) 2 0, 0, HI- Even if ini- 
tially this condition is not satisfied so that the system 
is in the topologically trivial state, the renormalization 
A — > A e g may make a topological phase possible since 
A c ff < A. Thus, periodic modulation of the chemical 
potential provides a way to tune the topological phase 
and so realize MFs by varying the parameter K/io. The 
rescaling Eq. holds only, of course, to the extent that 
off-diagonal couplings can be neglected; we address the 
generic case numerically below and show that more dra- 
matic changes in behavior are entirely possible. 

For periodic modulation of the Zeeman field, a sim- 
ilar analysis can be carried out by adding H z (t) = 
Kcos(ivt) X)i("it — to "0- Since only the Rashba 
term in Eq. ([3]) does not commute with the Zeeman term, 
the spin-orbit parameter is modified in the effective Flo- 
quet Hamiltonian: Ago — > ^SqJq(2K/u}). Thus, periodic 
Zeeman modulation cannot induce a topological phase 
transition if one keeps only the zero photon sector. How- 
ever, numerical investigation beyond the RWA [keeping 
all off-diagonal blocks of the effective Floquet Hamilto- 
nian cx J m - m '(2K/uj)} reveals that FMFs do, in fact, 
appear, and so we now turn to our numerical results. 
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FIG. 1: (color online) Quasi-energy spectrum for square-wave 
driven chemical potential. Parameters: r\ — 1.5 (full band- 
width D = 4rj = 6.0), V z = 1.0, A = 1.0, A so = 1-2, and 
(/Lti + /Lt2)/2 = 0.5. [Left panels]: quasi-energy near e = 0, as 
a function of driving period T for A/x = — /X2I = 2.6 (up- 
per), and as a function of driving amplitude A/j, for T = 1.75 
(lower). [Right upper panel]: quasi-energy near e = w/2 as a 
function of driving amplitude A/x for T = 0.6. [Right lower 
panel]: Finite size splitting (indicating the coupling between 
two FMFs at the two ends) for e — mode as a function of the 
number of sites in the chain (T = 1.75, Afx = 2.0). The finite 
size splitting shows exponential suppression accompanied by 
oscillations. There are N = 260 sites in the chain. Note: the 
unit used for the quasi-energies is uj/2 = n/T. 



For numerical convenience we consider square-wave 
driving of the chemical potential or Zeeman field: fi = \i\ 
for nT < t < (n + 1/2)T, and fj, = fi 2 for (n + 1/2)T < 
t < (n + 1)T (with n = 0,1,2,...), and similarly for 
V z . The evolution operator for the full period reads then 

_ ■ H 2 T _ .HjT_ 

U (T, 0) = e 1 2r e 4 , and the quasi-energy spectrum 
e a is obtained numerically using Eq. ([T]). In all cases here, 
the parameters at any instantaneous time correspond in 
the static system to the topologically trivial phase. 

The results for periodically modulated chemical po- 
tential are shown in Fig. Q] Clearly, one obtains stable 
e = Floquet Majorana zero modes (left panels) for a 
large range of parameters, as well as e = lo/2 modes (up- 
per right panel) [32[. Note that the parameters used in 
Fig. [T] are very far from those for which the RWA re- 
sult Eq. © yields a FMF: here V* - (p + 2r?) 2 < at 
all times, so no renormalized A can yield a non-trivial 
phase. Nevertheless, FMF appear once Afj, surpasses a 
threshold Afj, c . The figure shows that the threshold for 
an e = uj/2 FMF can be very small compared to that for 
an 6 = FMF, and also that the quasi-energy gap can be 
tuned by varying A/i. The splitting of a e = mode due 
to finite size effects is plotted in the right lower panel; 
it shows the expected decay of the level splitting as the 
number of sites, and hence the separation between the 
two FMF, increases. 



0.1 r 

0.05 



<3 -0.05 



-0.1 1 







ij;|: 


MM! 


!!! 

li! 


llin-jlil! 


iiili 


iii 
I 




ll'llllll 



0.1 



0.05 



-0.05- 



1.5 2 
Driving Period : T 



0.5 1 
Driving Amplitude 



1.5 
AV 



FIG. 2: (color online) Quasi-energy spectrum for square-wave 
driving of the Zeeman splitting, near e = 0. Parameters: 
rj = 1.4 (full band-width D = 4rj = 5.6), V 2 = 1.0, A = 2.0, 
Aso = 1-5, and (/zi + /tta)/2 = 0.0. Left panel: quasi-energy 
as a function of driving period T, for AV Z — \V z i — V^a| = 1-8. 
Right panel: quasi-energy as a function of driving amplitude 
AV Z , for T = 1.1. There are 260 sites in the chain. 

The quasi-energy spectrum with periodic Zeeman 
splitting is shown in Fig. [2] It also reveals FMFs. We 
stress once again that to obtain FMF in this case, the 
RWA is not enough and off-diagonal blocks of the Flo- 
quet Hamiltonian are crucial. 

Floquet Topological Qubit and Non-Abelian Statistics 
— A natural question is whether FMFs can form topo- 
logical qubits, as their static counterparts do. FMF can 
certainly encode quantum information: an FMF exists 
at all instantaneous times, and neither chemical poten- 
tial driving nor Zeeman driving changes the total fermion 
parity. Then, a more difficult question is whether FMFs 
obey non-Abelian statistics. We will provide an argu- 
ment for a 2D system, which can then be generalized to 
a ID network following the argument for static MF [29| . 

Suppose that FMFs are moved (which can be achieved 
by tuning the driving potential on and off, or chang- 
ing the driving amplitude,) along a path R(t) with the 
Schrodinger equation [H(R(t),t) - id t ]\$(t)) = 0. The 
position of the FMF R(t) is assumed to vary on a very 
slow time scale compared to the fast periodic driving. 
Then, it is convenient to separate the fast and slow time 
scales, and apply the two-time formalism of Floquet the- 
ory 0, §]: idt —> idt + id T , where t indicates the fast 
time and r denotes the slow time. Then the Schrodinger 
equation becomes 

id T \$(R(T), t)) = [H(R(r),t) - id t ] MR(t), t)>. (8) 

When we consider the dynamics on the slow time scale, 
the fast time t can be considered as aparameter. It was 
pointed out by Brcucr and Holthaus [7J (also see a recent 
discussion [30]) that a Floquet system follows a general- 
ized adiabatic theorem. Define the instantaneous (for r) 
quasi-energy eigenstates using the Floquet operator 



H{R(r),t)-id t \ \4> a (R( T ),t)) =e a (R(T))\(f> a (R(T),t)). 

(9) 

Suppose the system is initially in a Floquet state 
|$(i?(r = 0),t)) = \(f> a {R{T = 0),*)). Standard proce- 
dures in quantum mechanics can be applied to Floquet 



4 



states as long as the extended inner product mentioned 
above, ((•!•)), is used. Second order perturbation theory 
then yields 0, 

mWr),t)) = e-W.We-'x.Wfl^^r),*)) 

- \MR(r), t)) ««^»$& m ) , do) 

where 9 a (r) = dr'e Q ,(i?(r')) is the dynamical phase, 
and X a(r) = £ dT J {{<f> a {R(/))\id T >\<l> a (R(7 J )))) is the 
generalized Berry phase. Therefore, to avoid transitions 
to other quasi-energy states, the change in time scale r 
must be slow and the quasi-energy gap should be large: 
Mi?(r)) - e a (R(r))\ » |((^(i2(r))|0 T |^(J?(r))))|. 
We assume this condition is satisfied so that the system 
will stay in its initial Floquet state. 

The Floquet Majorana excitations can be described by 
a Bogoliubov quasi-particlc operator, 

7+(t)= /dr[ U (r,i?(T),0^ t (r)+ V (r, J R(r),^(r)], (11) 



where "i/j ( r ) (V>(r)) creates (annihilates) a fermion at r, 
and v = u* for a MF. A U (1) gauge transformation which 
changes the superconducting order parameter phase by 
2-7T [3j is allowed by using the extended space of the Flo- 
quet system [32[ . This causes a minus sign on both $ (r) 
and V>(r), changing the sign of the FMF operator as well. 
Due to such multivaluedness, a branch cut is necessary to 
define the phase of the wave function. So, the exchange 
of two FMFs 7i(£) and Jj(t) can induce a transforma- 
tion: %(t) — > 7j(i) and jj(t) — 7(f) (since one of 
the FMF, say 7/(f), must pass through the branch cut). 
For a ID network, the exchange of two FMFs (through 
a T-junction, for instance) flips the sign of the supercon- 
ducting pairing term, which results in exactly the same 
transformation as in the 2D case [2!| . 

Given two FMFs 71(f) and 72(f), one can form a 
non-local regular fermion d){t) = (71(f) + ij 2 (t))/^/2. 
Let |G(f)) be the Floquet BCS state which is annihi- 
lated by any Floquet quasi-energy operators. \G(t)) and 
d'(t)\G(t)} form a two-fold degenerate space. The ex- 
change of two MFs results in \G(t)) -> e iv \G(t)) and 
${t)\G(t)) e l v e ™l 2 d\t)\G{t)). The tt/2 phase differ- 
ence after the transformation signifies non-Abclian statis- 
tics [H 0]. 

The exchange of two MF can also induce an extra uni- 
tary evolution involving a 11011- Abelian Berry matrix Q. 
The form of the matrix can be generalized to a Floquet 
system [32| by replacing (-|-) with ((-|-)); the unitary evo- 
lution then reads 



U(t) = Pexp 



1 1 M(r')dr' 




(12) 



where P denotes path-ordering and M a| g(r) = 
((<j>a(R(T))\id T \cj>p(R(T)))) is the generalized non- 
Abelian Berry matrix |32j. We want to test whether 



M Qi a causes any extra phase difference that breaks the 
non-Abelian statistics of FMFs. First, the non-diagonal 
matrix elements of M. a p are zero since fermion parity is 
conserved (as emphasized above this is true for all driving 
scenarios). Second, we follow a procedure similar to that 
for a stationary MF [33 , 34 1 where the odd parity element 
i((G\d d T [d)\G)Y) is written as the sum of the even parity 
element i((G\d T \G)) and an extra term i((G\(dd T d<)\G)) . 
It is just this term that might affect the the phase dif- 
ference 7r/2 and so the non-Abelian statistics. By using 
Eq. (fTTj) and the MF condition Vi — u* one finds 



((G\(dd T S)\G)) = 



dt c?rRe(u 1 9 T U2 — u^tUi). 



(13) 

This term is exponentially small since it contains an over- 
lap of wave functions for spatially separated MFs; fur- 
thermore, it actually vanishes in the adiabatic limit. We 
conclude that the non-Abelian Berry phase does not de- 
stroy the desired statistics of FMFs. 

Summary — Periodic modulation of the chemical po- 
tential or the Zeeman field appears to be a promising 
way to produce FMFs, both of which can be realized in 
ID cold atom condensates. We find that Floquet MFs 
are robust and can be generated in a wide parameter 
range. This system may have potential for topological 
quantum computation since FMFs obey the same non- 
Abelian statistics as their equilibrium counterparts. 

D.E.L and A.L. acknowledge support from Michigan 
State University. The work at Duke was supported by 
US DOE, Division of Materials Sciences and Engineering, 
under Grant No. lde-sc0005237l 



[1] 
[2] 

[4] 

[5] 
[6] 

l<] 

[8] 

[9] 

[10] 

[11] 
[12] 



Electronic address: delphy@msu.edu 
G. Moore and N. Read, Nuc. Phys. B 360, 362 (1991). 
N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 
L. Fu and C. L. Kane, Phys. Rev. Lett. 100, 096407 
(2008). 

J. D. Sau, R. M. Lutchyn, S. Tewari, and S. Das Sarma, 
Phys. Rev. Lett. 104, 040502 (2010). 
J. Alicea, Phys. Rev. B 81, 125318 (2010). 
R. M. Lutchyn, J. D. Sau, and S. Das Sarma, Phys. Rev. 
Lett. 105, 077001 (2010). 

Y. Oreg, G. Refael, and F. von Oppen, Phys. Rev. Lett. 
105, 177002 (2010). 

M. Sato, Y. Takahashi, and S. Fujimoto, Phys. Rev. Lett. 
103, 020401 (2009). 

S.-L. Zhu, L.-B. Shao, Z. D. Wang, and L.-M. Duan, Phys. 
Rev. Lett. 106, 100404 (2011). 

L. Jiang, T. Kitagawa, J. Alicea, A. R. Akhmerov, D. 
Pekker, G. Refael, J. I. Cirac, E. Demler, M. D. Lukin, 
and P. Zoller, Phys. Rev. Lett. 106, 220402 (2011). 
K. T. Law, P. A. Lee, and T. K. Ng, Phys. Rev. Lett. 
103, 237001 (2009). 

J. D. Sau, S. Tewari, R. M. Lutchyn, T. D. Stanescu, and 
S. Das Sarma, Phys. Rev. B 82, 214509 (2010). 



5 



[13] M. Wimmer, A. R. Akhmerov, J. P. Dahlhaus and C. W. 

J. Beenakker, New J. Phys. 13, 053016 (2011). 
[14] D. E. Liu and H. U. Baranger, Phys. Rev. B 84, 201308 

(2011) . 

[15] A. Y. Kitaev, Phys. Usp. 44, 131 (2001). 
[16] V. Mourik, K. Zuo, S. M. Frolov, S. R. Plissard, E. P. A. 
M. Bakkers, and L. P. Kouwenhoven, Science 336, 1003 

(2012) . 

[17] J. R. Williams, A. J. Bestwick, P. Gallagher, S. S. Hong, 

Y. Cui, A. S. Bleich, J. G. Analytis, I. R. Fisher, D. 

Goldhaber-Gordon, Phys. Rev. Lett. 109, 056803 (2012). 
[18] L. P. Rokhinson, X. Liu, and J. K. Furdyna, 

I arXiv: 1204.42121 (2012). 
[19] J.-i. Inoue and A. Tanaka, Phys. Rev. Lett. 105, 017401 

(2010). 

[20] T. Kitagawa, E. Berg, M. Rudner, and E. Dernier, Phys. 

Rev. B 82, 235114 (2010). 
[21] N. H. Lindner, G. Refael, and V. Galitski, Nat. Phys. 7, 

490 (2011). 

[22] H. Breuer and M. Holthaus, Phys. Lett. A 140, 507 
(1989). 



[23] S. C. Althorpea, D. J. Kouria, D. K. Hoffmanb, and N. 

Moiseyev, Chemical Physics 217, 289 (1997). 
[24] J. H. Shirley, Phys. Rev. 138, B979 (1965). 
[25] The chemical potential in Ref. @, is /i — > /x + 2?7. 
[26] A. Eckardt, C. Weiss, and M. Holthaus, Phys. Rev. Lett. 

95, 260404 (2005). 

[27] H. Sambe, Phys. Rev. A 7, 2203 (1973). 

[28] A. A. Reynoso and D. Frustaglia, larXiv:1208.2742V 2 

(2012). 

[29] J. Alicea, Y. Oreg, G. Refael, F. von Oppen, andM. P. 

A. Fisher, Nat. Phys. 7, 412 (2011). 
[30] T. Kitagawa, T. Oka, A. Brataas, L. Fu, and E. Demler, 

Phys. Rev. B 84, 235108 (2011). 
[31] D. A. Ivanov, Phys. Rev. Lett. 86, 268 (2001). 
[32] See supplementary materials for further details. 
[33] M. Stone and S.-B. Chung, Phys. Rev. B 73, 014505 

(2006). 

[34] M. Cheng, V. Galitski, and S. Das Sarma, Phys. Rev. B 

84, 104529 (2011). 
[35] F. Wilczek and A. Zee, Phys. Rev. Lett. 52, 2111 (1984). 



I 

SUPPLEMENTARY INFORMATION 



In this supplementary information we (i) provide more details about the rotating frame analysis (specifically, the 
derivation of Eq. (7) in the main text), (ii) provide more data on the quasi-energy spectrum, (hi) discuss U(l) gauge 
invariancc in the extended space of Floquet system, and (iv) develop a generalization of the non-Abelian Berry matrix 
to Floquet system. 



Derivation of Eq. (7) in the main text 

We want to calculate matrix elements of the effective Floquet Hamiltonian in a basis of the rotating frame [l[ . The 
starting point is Eq. (6) in the main text: 



#Fioquet = m\H(t) - id t |{n<}; m')) 



(nn+na) {H + mui) e - iKsi :^ E i ("n+"ii)|{ n ^}) e -'(™-'-"'M (314) 



where Hq is shown in Eq. (3) in the main text. We note that the operator 5Zi(^»T + ^4-) f aus t° commute only with 
the superconducting term 



cj, + h.c. 



(S15) 



Therefore, we have 

#Fioquct = Smm> [({ni}\H (A = 0)|{n-}) + mu] 
1 r T 



S mm > [<{M|#o(A = 0)|K» + mu] 



1 

'TJo 



dt({m}\H< 



sc 



2K sin(wt) \ f AH 1 / . 2K sin(wi) 



so 



2! 



SC 



lK»e 



-i(m— m')u)t 



5mm' [({ ni }\H {A = 0)|K}) + mu] 



+- I dt({ni}\(Hsccos 



2K . , : 

sm(cjt) 



iHic sin 



. , s 

sm(cjf) 
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IK}) 



f\ \ — i(m— m )u)t 



(S16) 
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FIG. S3: (color online) Quasi-energy spectrum ear e = for square-wave driven chemical potential as a function of driving 
period T. Parameters: r\ = 1.5 (full band-width D = 4rj = 6.0), V z = 1.0, A = 1.0, A S o = 1.2, and (/n + fi 2 )/2 = 0.5. Upper 
panel: A[i — \fii — fx^\ — 1.4 (left), A/i = 2.0 (right). Lower panel: A^i = 2.6 (left), Afi — 3.2 (right). There are 260 sites in 
the chain. Note: the energy unit for the quasi-energies is cj/2 = n/T. 

where H^j? is an anti-Hermitian operator 

H£? = AY,(4A- h - c )- ( S17 ) 

i 

Since different photon sectors (labeled by m) are separated by the energy u>, then in the limit uj D (D is the 
band- width) , we need consider only the zero photon sector, m = 0. This is the rotating wave approximation (RWA). 
Then, the effective Floquet Hamiltonian becomes 

({m}\H Q (A = 0) + J (2K/u})H sc \{n' i }) (S18) 

Here, Jo(x) denotes zero order Bessel Function of the first type. The effective Floquet Hamiltonian within RWA 
has exactly the same form as Hq in Eq. (3) in the main text except that the pairing coupling A is renormalized to 
A c ff = A Jq(2K/ui). Note that the integral of the second term, iHgc , is zero for m = m! . As we emphasized in the 
text, non-diagonal terms, m 7^ m', proportional to J m -m'i may be important. We have given one example for the 
case of Zeeman field modulation. 



More data for Quasi-energy spectrum 

In the main text we showed a few representative examples of the quasi-energy spectrum as a function of driving 
period or amplitude. In this section we elaborate on this analysis and uncover a rather complicated structure of the 
quasi-energy states. This spectrum was found by numerical diagonalization keeping all the non-diagonal terms in the 
effective Floquet Hamiltonian. 

Fig |S3l shows the quasi-energy spectrum as a function of driving period T for square-wave chemical potential 
driving. As the driving amplitude increases from A/u = 1.4 to 2.0, the Floquet MF in the region T E (1.4, 1.9) appears 
gradually, and the quasi-energy gap becomes larger for the FMF in the region T E (1.0, 1.35). However, as the driving 
amplitude further increases, the quasi-energy gap becomes smaller (from A/z = 2.0 to 2.6) for FMFs in the region 
T € (1.0, 1.35). The T g (1.0, 1.35) FMFs are even be killed for A^i = 3.2. 

Fig. IS4I shows the the quasi-energy spectrum near e = lu/2 with periodic chemical potential as a function of driving 
period T. Clearly, one can see the e = uj/2 FMF mode appears in large region of parameter space. The quasi-energy 
spectrum near e = uj/2 with periodic Zeeman splitting is shown in Fig. IS51 which also reveals FMFs. 

Fig. IS6I shows finite size splitting, which indicates the coupling between two FMFs at wire ends, as a function of 
1/N for the e = mode. N is the number of sites in the chain. The splitting of two FMFs should decay exponentially 



7 



as the length of the wire increase: ~ exp(— L/£), where the length of the wire L = aN with lattice spacing a, and 
£ is the superconducting coherent length. We fit the data using the function form a ■ exp(—N/b) + c. As N — s- oo, 
the finite size splitting goes to zero for A/i = 1.9, 1.8, and goes to finite values for Apt = 1.7, 1.6, which shows the 
threshold A/j. c for the FMF is larger than 1.7 and smaller than 1.8. 

One important quantity in Floquet theory is the mean energy Q corresponding to the average expectation value of 
Hamiltonian over a full period: E a = (l/T) dt(<p a (t)\H(t)\4> a (t)) = e a — ujde a /duj. Numerical data shown indicate 
that the partial derivative part of this expression vanishes for MF zero modes (except near the transition point), so 
Eq = eo = and = = w/2. Once weak heating effects are considered, one expects that the driven system 
tends to occupy the Floquet state with the lowest mean energy [j| ||. Therefore, we focus on the e = FMF, that 
have the lowest mean energy. 



U(l) gauge invariance in the extended space of Floquet system 



For the system described by the BdG Hamiltonian, there is a U(l) gauge invariance Q: if the global phase of the 
superconducting order parameter A is shifted by <f>, i.e. A — > Ae 1 ^ , it is equivalent to rotating electron creation and 
annihilation operator ip — > e^^tp and ip^ — > e — ll ^/ 2 -0t_ 

Let us consider the instantaneous quasi-energy eigenstates (also shown in Eq. (9) in the main text) 



H{R{r),t)-id t \4> a [R{ T ),t)) =e a (R(T))\<f> a {R{T),t)). 



(S19) 



The quasi-energy state of the Floquet system is defined in a extended Hilbert space Q , therefore, one can expand the 
operator H(R(r),t) and state |0 a (i?(r), t)} as the sum of different photon-sector Q 



H(R(r),t) - ^ e" 

n 

\MR(r),t)) - 5> 



f ^(i?(r)), 
V«n(^(r))). 



(S20) 
(S21) 



Then, the Floquet Hamiltonian of the extended Hilbert space can be written as 

/■•::: \ 



H 



Hq + LU H± H2 

H-i H Hi 
H—2 H-i Hq — uj 



V 
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It is easy to check that the £7(1) invariance exists for all the matrix elements: • • • -H_2, H—i, Hq, H\,H2, 
the U(l) gauge invariance also exists for the Floquet system defined in the extended Hilbert space. 



(S22) 



Therefore, 
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FIG. S4: (color online) Quasi-energy spectrum near e = cj/2 for square-wave driving chemical potential. Parameters: r\ — 1.5 
(full band-width D = 4rj — 6.0), V z — 1.0, A = 1.0, Aso = 1-2, and (^i + /i2)/2 = 0.5. The spectrum is as a function of driving 
period T for A/i = - fi 2 \ = 2.0 (left), and for A/i = 2.6 (right). N = 260. 
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FIG. S5: (color online) Quasi-energy spectrum near e = ui/2 for square- wave driving of the Zeeman splitting. Parameters: 
r\ = 1.4 (full band-width D = 4rj = 5.6), V z = 1.0, A = 2.0, A SO = 1.5, and (m + fi 2 )/2 = 0.0. N = 260. Left panel: 
quasi-energy as a function of driving period T for A V z = 1.4. Right panel: quasi-energy as a function of driving amplitude 
AV Z for T = 1.25. 




FIG. S6: (color online) Finite size splitting (indicating the coupling between two FMFs of two ends) for e = mode as a 
function of 1/N, where indicates the number of sites in the chain. Parameters: rj — 1.5 (full band-width D — 4rj — 6.0), 
V z = 1.0, A = 1.0, Aso = 1.2, and (fii + fJ,2)/2 = 0.5. The period is T = 1.75. The curves from bottom to top are for 
Afi — 1.9, 1.8, 1.7, 1.6. The lines show the best fit to the function a ■ exp(—N/b) + c. 

The Floquet MF excitations shown in Eq.(10) in the main text can also be written as 

7iW = E e " i "^n, (S23) 

n 

where 

7L = I dr[u n (r,R(T))ft(r)+v n (r,R(T)Mr)]. (S24) 

When the over all phase of the superconducting order parameter change 2tt, the Floquet MF excitation change sign: 
iL -> -itn and thus jl(t) -> -7* (i). 



Non-Abelian Berry Matrix for Floquet System 



As shown in Eq.(8) in the main text, the adiabatic evolution of the Floquet system can be described by a two-time 
Sclirodinger equation 0,il (also see Eq.(8) in the main text) 



id T \$(R(T),t)) = H(R(r),t) - id t \$(R(r),t)), 



(S25) 
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Here, we will show if the quasi-encrgy degeneracy occurs in Floquet system, the adiabatic evolution of the Floquet 
system can be described by a generalized non-Abelian Berry matrix, as in the static counterpart Q. 
Consider an instantaneous quasi-energy equation for a fc-fold degenerate Floquet states 

H(R(T),t)\(t> a (R(r),t)) = e(T)\^ a (R(T),t)) , (S26) 

where a = 1,2, ■ • • k. If the quasi-energy difference is large between this subspace and other states, the transitions to 
the states outside the A:-fold subspace can be neglected within the adiabatic approximation. For any time r, the wave 
function of the system can be decomposed into a linear combination of the Floquet quasi-energy state 

MR(r),t)) =e- i f°^W£c a (T)\<f> a (R(r),t)) (S27) 

a 

Feed Eq. (|S27[) into Eq. (|S25|) . one can obtain 

^(i9 r c Q (r))|0 Q (i?(r),t)) +J2c a (T)i8 T \MR(r),t)) = 0, (S28) 

a a 

Projecting the equation to the state ((f>p(R(T),t)\ and carrying out the integral (1/T) J Q T dt, one finds 

id T cp{T) = -J2 M Pa c a (T) (S29) 

a 

where 

i - T 



M Pa { T ) = - dt{<Pp{R{r),t)\id T \<i> a {R{T),t)) (S30) 
1 Jo 

corresponds to elements of the fc-by-fc matrix M(r). Then, it is easy to check if the system is initially in the Floquet 
state \c/> a (R(T = 0), t)), the time (t) evolution of such state can be written as 

\* a (R(T),t)) = U(r)|0 Q (i?(r = 0),t)>, (S31) 

where 



U(t) = Pexp 



M(r')dr' 



(S32) 



and P denotes the path-ordering. This is the evolution operator given in Eq.(12) of the main text. 
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